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ABSTRACT 

The measurements of the Doppler shifts of the Fraunhofer hnes, scattered by the dust 
grains in the solar F-corona, provides the insight on the velocity field of the dust and 
hence on its origin. We report on such measurements obtained during the total eclipse 
of March 29, 2006. We used a Fabry-Perot interferometer with the FOV of 5.9° and 
the spectral resolution of about 5000 to record Fraunhofer spectral lines scattered by 
the dust of the F-Corona. The spectral region was centered on the Mgi 5172. 69A line. 
The measured line-on-sight velocities with the amplitude in the range from -10 to 
10 km-s~^ show that during our observations the dust grains were on the orbit with a 
retrograde motion in a plane nearly perpendicular to the ecliptics. This indicates their 
cometary origin. Indeed, at the end of March, 2006, SOHO recorded several sungrazing 
comets with the orbital elements close to what was deduced from our measurements. 
We conclude that the contribution of comets to the dust content in the region close to 
the Sun can be more important albeit variable in time. We also deduce that the size 
of the most of the dust grains during our observations was less than 0.1 |J.m. 
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1 INTRODUCTION 



The work of [van es tablished, after an earlier 

suggestion by Grotrian (as cited bv lvan de Hulst|[T947l '). that 
the emission of the solar outer corona, called F-corona after 
the Fraunhofer lines composing its spectrum, is due to the 
scattering of the solar hght by the dust particles. 

In recent decennia it became more clear that the dust 
component near the Sun cannot be described as a mere ex- 
trapolation of the zodiacal light disc with the dust grains 
orbiting in the plane of the ecliptics and slowly drifting to 
the Sun under the Poynting-Robertson effect as one could 
accept in a first approximation. 

The measurements of the Doppler shifts of the Fraun- 
hofer lines in the spectrum of the F- corona, first done durin g 



the total eclipse of July 31, 1981 by Shche glov et all (1 19871), 
and t he analysis of the derived velocities by Shestakoval 
1 19871 ). indicated that, although most of the dust particles 
are orbiting in the ecliptic plane, there are also grains orbit- 
ing at high values of inclination angle i, implying that the 
region close to the Sun receives a dust contribution from 
long-period comets. Also, the in situ measurements of the 
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dust particle velocity distribution on board of HELIOS-1 
showed the existence of two distinct comp onents in the inter - 
planetary dust, cometary and asteroidal l|Griin et al.lll980l ). 
Among more recent results, the high frequ ency of the Sun- 
grazi ng comets recorded by SOHO (M acQueen fc St. CvJ 
Il99ll ) makes also to consider the cometary contribution 
as important. Furthermore, studies of extrasolar planetary 
and dust systems showed the exist ence of Falling Eva porat- 
ing Bodies, possibly comets (e.g. iBeust et al]|l998l ). This 
observational progress was followed by a thorough theo- 
retical modelling taking into account co mplex composition 
of the dust in the s olar vicinity (see iMann et al.l |2000| . 
iKobavashi et al.ll2008l ). 

It would have been important to obtain further mea- 
surements of the line-on-sight (thereafter LOS) velocities of 
the dust near the Sun to have a better idea of the rela- 
tive contributions of the cometary and the ecliptic disc dust 
components and of their possible variation in time. After the 
first work of 1981, these measurements, to the best of our 
knowledge, were successfull y attempted only onc e, during 
the eclipse of July 11, 1991 (|Aimanov et al.l 19951 ) confirm- 
ing the first results. 
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Here we present new higher quality measurements of the 
LOS velocity of the dust in the F-corona obtained during the 
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Figure 1. The Fabry-Perot interferograms of the dayhght sky and of the circumsolar region up to elongation e ~ 11. Here, both axes 
are in pixels, the platescale is 20.66 arcsec/pix, the FOV= 5.9°. From left to right: a comparison daylight sky frame "D", the total 
eclipse frame "E" , and the frame "L" taken at the end of the totality, but already including the green coronal emission line of [Fe XIV] 
at A5302.86A. The fiducial wire (black line) indicates the East- West direction, the West is to the left along the wire and the North is 
down. The frames are rotated so that the horizontal axis of the frames coincides with the ecliptic plane. The frames are dark current 
subtracted. For more details see the text. 



eclipse on March 29, 2006 fsee ^Shestakov a et al ] |2007l . for a 
preliminary report), and discuss their implications. 



2 OBSERVATIONS 

The site of the observations was the village of Mugalzhar in 
the Aktobe region of the Republic of Kazakhstan, ai (f} — 
48° 35' and A = 58° 27', situated in the middle of the totality 
band. According to our calculations, the beginning of the 
total eclipse was at UT llh 32m 40, its end at llh 35m 30s, 
and the duration of the total phase was 170 s. The Sun was 
at 27? 5 above the horizon. The weather was mostly rainy 
and windy during the days preceding the eclipse, but on the 
eclipse day, the sky was clear with no wind and excellent 
transparency. 

We used the s ame optical set-up as in our previous mea- 
surements in 1991 (IShcheglov et al.lll987^ ■ i.e. a Fabry-Perot 
(hereafter FP) spectrometer with a coronographic mask re- 
jecting the hght of the solar corona and thus reducing the 
background. The entrance lens of 10 cm diameter has its 
focal plane at the field lense, the latter is followed by a 
collimator. The FP etalon in series with an order separa- 
tion filter is placed in the parallel beam close the exit pupil 
imaged by the field lense. It is followed by a photographic 
objective and a CCD detector. The latter was an Apogee 
Alta-10 CCD device with 2048^ pixels of 14 (j.m size. During 
the data reduction, the frames were binned by 2 x 2 pixels, 
so that to avoid any confusion we will refer from now on 
to the detector format of 1024^ pixels with the 28 \i.m pixel 
size. 

The solar angular radius at the time of the eclipse was 
7^0 = 961 arcsec, the measured on the CCD value was 
46.5 pix. This gives the equivalent focal length of 279.4 mm, 
and the platescale of 20.66 arcsec/pixel. The field of view 
(FOV) is 5.9°, corresponding to the elongations e < 11, 
however, to reduce the background hght, the central region 



around the Sun, e < 2.6, is hidden by the coronographic 
mask (hereafter, the elongation e is given in angular solar 
radii 7^.0 ). 



2.1 Spectro-imaging 

We used a mechanically adjustable FP etalon with the free 
space A = 70 vim. The order separation filter has the FWHM 
= 10 A with the central wavelength close to that of the Mg I 
line at 5172.69 A. In the wavelength space, the set-up is such 
that a frame includes 4 FP orders. The highest, on-the-axis 
order = A/A, is at the centre of the frame and is unseen 
due to the use of the coronographic mask; the orders lower 
than — 4 are outside of the FOV. 

The dark subtracted frames obtained during the eclipse 
are shown in Fig.[T] The left frame, denoted as "D" , is one of 
the calibration interferograms of the daylight sky, scattered 
on a white screen, recorded shortly before and after the total 
phase of the eclipse. The use of the white screen allows a 
homogenous brightness distribution over the field of view 
without changing the spectrometer position which remains 
at the same position as during the total phase. The frame 
shows concentric rings of the Fraunhofer absorption lines 
scattered in the terrestrial atmosphere. 

The dark circle at the centre of the frame corresponds 
to the coronographic mask. The dark line corresponds to a 
fiducial wire indicating the East- West axis, the West to the 
left along the wire and the North down. For convenience 
of interpretation, the frames are rotated so that the frames 
horizontal axis coincides with the ecliptic plane. The frames 
are dark current subtracted. 

The frame "E"in the middle of the Fig. [1] is our main 
scientific exposure, of 130 s, taken on the circumsolar region 
during the eclipse. It was started a few seconds after the 
beginning of the totality and stopped well before its end, so 
that contamination by coronal or chromospheric lines was 
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Figure 2. The interferogram of the circumsolar region during the total eclipse in different states of data reduction procedure. Left: after 
median filtering and dividing on the flat-field as described in the text. Right: The same with the low frequency two-dimensional trend 
subtracted. The elongation e is in angular solar radii 7?.q . 



avoided. The concentric dark rings are the solar Fraunhofer 
hnes scattered on the dust particles of the F-corona. 

Finally, the frame "L"at the right of the Fig.[T] with 
the exposure time of 20 s was started close to the end of 
the total phase and lasted a few seconds beyond, so that it 
caught the Ught of the solar corona which was just appear- 
ing from behind the mask (but not yet the photosphere). 
On this frame, one can see, additionally to the absorption 
spectral features of the F-corona, the bright emission rings 
of the green coronal [Fe XIV] line at A5302.86A. Its wave- 
length lies far from the centre of the our narrow-band filter, 
however the emission is so strong that its light passed in 
the filter transmission wings. The corresponding interfer- 
ence rings well visible in Fig. [1] traced the location of the 3 
used Fabry-Perot orders and gave useful reference points for 
the data reduction. Its measured FWHM = 1.2±0.lA gives 
the spectral resolution of the instrument. 

The daylight sky interferograms "D" provided the wave- 
length standard and allowed to eventually measure the 
Doppler shift of the dust particles of the F-corona, while 
the interferogram of the coronal green line A5302.86A al- 
lowed to calibrate the spectral geometry of the frames, and 
in particular to accurately define the rings centre as it will 
be discussed in Section[3l 

2.2 Brightness calibration 

The weather conditions being excellent and stable, the 
Fabry-Perot data were calibrated in brightness relative to 
that of the solar disc, Bq. Two 0.1s exposures on the Sun 
through a combination of a neutral and green filters were 
taken to the East, at e = 4.0, and to the West, at e = 4.4 
from the Sun position at the total phase of the eclipse. The 
Sun elevation for the calibration and for the eclipse frames 
was nearly the same, about 30° . The uncovered surface of the 



solar disc was 0.94 and 0.25 respectively for the two expo- 
sures. Taking this into account, the daylight sky brightness 
off the eclipse is estimated to has been 4 ■ 10~^-Bq. 

The brightness of the F-corona continuum emission. 
Bp, was measured to be 4.1 • IO'^Bq at the elongation 
e = 4.0, it was 3.9- 10~^Bq at e = 4.4, and 1.7- 10"^Bq for 
the e range from 9 to 10. The brightness decrement in the 
F -cor ona is in agreement with the results of lKoutchmv et al] 
(jl978i ) . Compared to the daylight sky, the F-corona is fainter 
by a factor of 10*. 



3 DATA REDUCTION 

The extraction of the Doppler shift of the solar absorption 
lines out of the Fabry-Perot interferograms needs a thorough 
data reduction process. It is worthwhile to be describe here 
in details. 



3.1 Noise and trends handling 

After the bias and dark current subtraction, all the frames 
were 2x2 binned, which increased the signal-to-noise ratio 
without changing spectral or spatial resolution. The next 
step was defining all regions meaningless for further reduc- 
tion, which includes the central dark region corresponding 
to the coronographic mask, that of the fiducial wire, the part 
lying out of the field-of-view of the detector and the regions 
of strong light induced by spurious reflections. The resulting 
numerically masked area was about 30%. 

The "hot", "cold" and "dead" pixels were cleaned out 
using a median filter. The visual analysis showed that the 
"clouds" of defective pixels did not exceed a region of 10-12 
pixels a size. We used therefore a non-linear circular filter 





Figure 3. Extracted spectrograms: The flux integrated over the position angle PA in a function of the radius r. Upper plot: The daylight 
sky, frame "D" ; Middle plot: The circumsolar region during the total eclipse, frame "E" ; Lower plot: Same but also including the strong 
coronal emission [Fcxiv] A5302.86A, frame "L". The relevant Praunhofer lines and their wavelengths are indicated. Note that there are 
3 overlapping spectral orders, they are traced by the [Fexiv] emission pics. 



covering 31 neighbour pixels (i.e. 2n+l, where n is the meix- 
imum scale of "clouds" ) . The correction was applied only to 
the pixels with the count exceeding 20% of the median value. 
The number of corrected pixels was less than 4% on the 
eclipse frames and less than 0.5% on the calibration frames. 

For convenience of interpretation, all frames were then 
rotated so that the horizontal axis coincides with the ecliptic 
plane. 

The next step was the flat-fielding. Usually, it is done by 
a straightforward division by a averaged flat-field frame with 
a subsequent masking of the pixels resulting from the divi- 
sion by zero or close to zero values. We decided to optimize 



this operation by adding a small constant to the flat-field 
frame. This is similar in its spirit to the Tikhonov regulari- 
sation of ill-defined inverse problems. The value of the con- 
stant was defined by analysing the histogram of the counts 
of the flat-field frame. The aim of this operation is to keep 
moderate, or negligible, the change induced in the signal- 
to-noise ratio. Thus, the division on the flat-field does not 
change the structure of the frames and avoids a useless ad- 
dition of a noise. 

Further, we measured the two-dimensional low fre- 
quency trend of the resulting frames in order to subtract 
it and to keep only useful spectro-spatial variations. It was 
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done by using iteratively a linear circular moving average 
filter. The best filter diameter and the number of iterations 
were searched by trials in such a way that the frame with 
the subtracted 2D trend would still keep the structures on a 
10-15 pixels scale, which is that of the recorded Fraunhofer 
lines. The best filter had the diameter of 13 pixels, covering 
137 closest pixels, and the best number of iterations was 4. 
Higher the number of iterations, closer the iterative filter 
to a linear gaussian smoothing filter, exp[—{x^ + i/^/(2crg)], 
where Ug defines the degree of the smoothing. The advan- 
tage of using an iterative filter is the possibility to control the 
achieved smoothness by limiting the number of iterations. 

We tried square and rectangular moving average filters, 
however they give rise to spur line-like features. We also 
applied the methods of high and low enveloping curves, but 
it did not give better results. 

The FP interferogram after median filtering and sub- 
traction of the 2D low frequency trend is shown in Fig.[51 

Finally, the data were passed through the gaussian filter 



with G„ 



1.3 pixels, or FWHM ~ 3 pixels which is close to 



the measured FWHM of the point-spread function of the 
experiment. 



3.2 Defining the interferograms centre 

In an ideal Fabry-Perot interferogram, the wavelength A is 
constant on a circular ring, and varies with its radius as 
r^. Let Xa and Yc denote the coordinates of the centre of 
the interference rings. Their values depend on many optical 
parameters, which can vary with temperature and mechan- 
ical flexures, so that they must be carefully defined for each 
frame. After different trials, we found that the most reliable 
way to measure the centre position was to use a correlation 
method in the following way. 

Let us denote z = . We adopt a first guess of the 
centre coordinates Xq and Yq, and divide the frame on two 
sub-frames, left and right, symmetrically with the respect 
to Xq. For each sub- frame, we compute counts vs z, which 
gives us two functions L{z) and R{z) respectively for the left 
and the right sub-frames. We compute then the correlation 
of L and R , and vary Xq. The maximum of the correlation 
gives the value of Xq, which is adopted as X^. In a similar 
way, we define the best value of Yc, correlating the upper 
and the lower sub-frames. 

Such a correlation measurement can give false and bi- 
ased results if, for example, there is a strong asymmetry in 
the intensity of interference patterns. To have an additional 
check, we applied another method to smaller parts of the 
patterns, using only arcs of the rings. For a given ring, we 
take first guess values of Xc, Yc, r, and the width of the ring 
dr. Varying the values of Xc, Yc and r, we find the values 
such that the sum of counts is the less (for absorption lines). 
To insure a good statistics, the value of dr must be suffi- 
ciently large, but without covering neighbor ring patterns. 
We used values of dr in the range from 2 to 12 pixels. The 
difference of the centre coordinates defined by this method 
and that of correlations did not exceed 0.3 pix, which is at 
the level of the expected uncertainty. 

The resulting values of the center coordinates are de- 
fined with an accuracy better than 0.3 pix. 



3.3 Reduced spectrograms 

Once the centre of the interference rings was found, it is 
convenient to transform the data presentation from carte- 
sian to polar coordinates (r, cj>). The frames were oriented so 
that the polar angle <j) and the position angle on the sky, PA 
are the same for all of them. The extracted spectrograms in 
the form of the flux integrated over all values of the posi- 
tion angle PA in a function of the radius counted from the 
interferogram centre are plotted in Fig.[3] for the daylight 
sky, the circumsolar region during the total eclipse and the 
frame "L" including the coronal Fexiv, A5302.86A line. The 
3 pics of the coronal emission line trace the 3 Fabry-Perot 
spectral orders. The relevant scattered Fel, Fell and Mgl 
Fraunhofer lines, in absorption, are indicated. Their wave- 
lengths, the values of the ring radius r and the correspond- 
ing value of the elongation e are given in the Table[T] Albeit 
the used Fabry-Perot orders overlap, the spectral features, 
fortunately, are distinct and can be easily identified. 



3.4 The line-on-sight velocity of the dust. 

The Doppler shift between the eclipse and the daylight 
spectrograms was measured by cross-correlation. The cor- 
responding dust LOS velocity V writes: 



V = 



2/2 



or, substituting the numerical values: 

3 X 10*^ X 0.028^ X A? ^2 
^ - 2^86:65^ = ^ 



(1) 



(2) 



where / = 86.55 mm is the focal distance of the camera 



objective, 0.028 mm is the pixel size, = 



skyi 



where 



re and Vsky are the radii in pixels respectively for the eclipse 
and daylight sky frames. The pixels convert to the elongation 
e given that TZq = 46.5 pixels. 

The sampling in elongation e and in position angle PA 
was as follows. We used intervals of e as given in Tab.[2] 
choosing their spanning as regular as possible in a function 
of the useful spectral lines. The full range of PA was simply 
divided on 36 sectors 20° wide each. 

For convenience of the analysis, we computed two kinds 
of average velocities, K, which is the average over all values 
of PA, and Vpa, which is the average over all values of the 
elongation e. 

The value Vt is the function of e only; it is plotted in 
Fig.|3] If the LOS velocities of the dust grains are distributed 
in a central symmetry with the respect to the Sun, as it is 
expected for a pure circular motion, then Ve is 0; if not, it 
reflects the presence of a radial motion, an infall or outflow, 
at this particular value of e. The Fig.|4] shows that radial 
motions with velocities about 10 km ■ s~^ or less might be 
present at e = 4.5 and e « 9, and they are absent, or very 
small, in between these elongations. 

The value Vpa is a function of the position angle PA 
only, it is plotted in Fig.[Sl It is so particular that we give 
also for comparison a similar plot from the eclipse on July 
31, 1981. For 2006, the triangle marks indicate measure- 
ments using two different daylight sky interferograms; their 
scattering provides an estimate of the uncertainty 5VpA — 
1.6 km ■ s"\ The values for the PA range 360° - 450° are 



6 L. I. Shestakova et al. 



Table 1. Spectral lines identification. 



Radius 


Elongation e 


Wavelength 


Element 


(pixels) 


(7^0) 


A(A) 




182 


3.92 




TTp T 4- TT 


196 


4.22 


O ±U 1 .oo 


IVTd- t 

iVlgl 


220 


4.72 


5183.62 


Mgi 


238 


5.12 


5162.28 


Pel 


266 


5.72 


5159.06 


Fel 


298 


6.40 


5172.69 


Mgi 


328 


7.06 


5167.33 


Mgi 


344 


7.40 


5183.62 


Mgi 


358 


7.68 


5162.28 


Fei 


374 


8.04 


5159.06 


Fei 


402 


8.64 


5172.69 


Mgi 


425 


9.14 


5167.33 


Mgi 


448 


9.64 


5162.28 


Pel 


460 


9.88 


5159.06 


Pel 


485 


10.42 


5172.69 


Mgi 



Table 2. Areas used for the Doppler shift measurements. 



Elongation range 


Mean e 


Lines used 


3.12 - 4.09 


3.66 


5172.7, 5169.0 


3.12 - 4.82 


4.09 


5169.0, 5167.3 


4.09 - 4.82 


4.52 


5167.3, 5183.6 


4.82 - 5.63 


5.27 


5183.6, 5167.3 


5.63 - 6.37 


6.02 


5159.1, 5172.7 


6.36 - 7.01 


6.71 


5172.7, 5167.3 


7.01 - 7.74 


7.40 


5167.3, 5183.6, 5162.3 


7.74 - 8.43 


8.08 


5162.3, 5159.1 


8.43 - 9.12 


8.77 


5172.7, 5167.3 


9.12 - 9.63 


9.38 


5167.3, 5162.3 


9.63 - 10.25 


9.98 


5162.3, 5159.1, 5172.7 



added for convenience, they merely repeat those for — 90°. 
Plotted are also the least squares sinus fits in the form: 
VpA = fci + ■ sin{kz ■ PA + ^4). The fit coefficients, 
for 2006, are as follows: ki = 0.3 ± 0.3, = 18.7 ± 0.3, 
ka = 0.01773 ± 0.00016 and ki = 1.22 ± 0.04, and for 1981, 
fci = 8.8 ± 2.5, fc2 = 22.0 ± 3.5, ks = 0.019 ± 0.001 and 
k4 = 2.4 ±0.3. 



4 DISCUSSION 

For the prograde orbits being strictly in the ecliptic plane, 
the minimum LOS velocity Vpa should be at 90°, i.e. to the 
East from the Sun, and the maximum Vpa at 270°, i.e. to 
the West. The velocity curve measured on July 31, 1981, is 
close to what is expected from such a motion (with, possibly, 
a slight difference, the minimum of Vpa being at 120° and 
the maximum at 300°). 

As to the Vpa curve on March 29, 2006, the values of 
PA of the extrema differ dramatically from what is expected 
from the prograde motion in the ecliptics, namely the min- 
imum VpA is at 195°, and the maximum is at 15°, meaning 
a retrograde motion in a plane nearly perpendicular to the 
ecliptics. 

The dust rotating in the prograde direction in the eclip- 
tics is also barely present showing the "jump" of Vpa at 



PA ~ 120° and the symmetric to it another "jump" at 
PA ~ 300° (see Fig.El upper plot). 

Let us verify whether the the LOS velocity behavior 
agrees with what one would expect from the keplerian mo- 
tion. This can be seen from the variation of the LOS velocity 
with the elongation e at a given PA. We assume that the 
elongation e of the scattering dust and its heliocentric dis- 
tance are in a linear relation. We choose the values of 15° 
and 195° for the PA, which corresponds to the extrema of 
Vpa curve. First of all, the LOS velocities in the considered 
Pj4's were averaged within a sector of ±30° wide, then we 
subtracted from it the Ve value. The resulting "orbital" LOS 
velocity Vorb{^) is given in Fig.|4]together with the linear fit. 
The latter gives V'ort.(e) = —2.2 ± 3.1e with the la uncer- 
tainty of ±2 on the additive term and ±0.3 on the factor. 
The success of the fit indicates that the motion is indeed 
keplerian. 

Summarizing what is indicated by the Vpa curve on 
March 29, 2006, we conclude that the bulk of the dust grains 
was in a keplerian retrograde motion in a plane nearly per- 
pendicular to the ecliptics. In the inner Solar system, the 
only objects having this kind of motions are comets. 

It happens that indeed around the eclipse date, the 
SOHO spacecraft recor ded a series of sungrazin g comets 
dMazzucato et all I2OO6I) of th e Kreutz group ijMarsdenI 
1 19671 , ISekanina fc Chodasll2007l ). One of them, CK06F050, 
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Figure 4. The line-on-sight velocity of tlie dust Ve vs the elongation e; Ve is the average over all range of the position angle PA. 




Figure 5. The line-on-sight velocity of the dust averaged over the studied range of elongations, VpA, the position angle PA on March 
29, 2006 (upper plot) and on July 31, 1981 (lower plot) together with the fitted sinus curves. For 2006, the triangle marks indicate 
measurements using two different daylight sky interferograms. The values for the PA range 360° — 450° were added for convenience, 
they merely repeat those for — 90°. The PA values of the extrema of the fitted to Vd sinus curves are given. For the circular orbits 
lying strictly in the ecliptic plane, the extrema should be at 90° and 270° . 



fell onto the Sun a day before the eclipse, another one, 
CK06F070, fell a day and a half after the eclipse. For con- 
venience of the r eader , the orbital elements as given by 
iMazzucato et all l|2006l ). are reproduced in the Tab.[3l It is 
quite possible that there were also a chain of smaller frag- 
ments in between the entities recorded by SOHO, and thus 



the FOV of our eclipse frame is filled by the dust brought by 
the Kreutz group comets moving to the Sun on retrograde 
orbits. 

The date of our observations is 8 days past the equinox, 
so that the line-of-sight to the Sun projects to approximately 
8° East from that to the vernal point fi. The longitude of 
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Figure 6. Upper plot: The difference Vo^i, = V — Ve vs the elongation e in the directions of PA = 15° (rectangles) and PA = 195° 
(triangles) and the corresponding linear fits. The values of V are the average over ±30° wide sectors, they are given in absolute values. 
The actual sign is positive in the first direction and negative in the second one. Lower plot: The average of the two and the linear fit to 
it, V„rt = -2.2 + 3. le 



Table 3. Orbital elements of sungrazing comets at the end of March, 2006 



Comet name 


Epoch 


Perihelion q 


Perihelion 


Longitude of 


Inclination i 






in AU 


argument 


the node 




CK06F050 


March 21.96 


0.0050 


82.38 


4.20 


144.57 


CK06F060 


March 23.04 


0.0333 


56.09 


75.03 


74.13 


CK06F070 


March 28.64 


0.0050 


84.83 


3.42 


145.72 


CK06F080 


March 31.10 


0.0052 


84.02 


5.67 


144.58 



the node of the indicated comets is distant from the point, 
according to Tab.[3l by 3° — 6°. This means that our line- 
of-sight is inclined with the respect to the line of nodes by 
only 2° — 5°. The perihelion argument lies at 84° — 85° from 
the line of nodes. 

This means that our image plane is nearly perpendicular 
to the orbital plane of the comets, and the line-of-sight slides 
over it. 

If the dust grains were strictly confined to the orbital 
plane, we would see the selected direction indicated by the 
extrema of the LOS velocities at PA ~ 55°, which would 
correspond to the inclination of the parent comet's orbit 
of i = 145°. But our data indicate a different axis for 
the extrema velocities, namely PA ~ 15°, which indicates 
i = 105° , i.e. the orbital plane of the dust grains is turned 
with respect to that of the comets by 40° . The orbital plane 



cannot change under the gravitational force. Hence, we have 
to assume another reasons, the effect of the magnetic field 
on electrically charged dust particles being the most plausi- 
ble. It follows that the size of particles is quite small, which, 
in turn, suggests their cometary origin. According to the 
det ailed models of the dust grains dynamics near the Sun 
by (|Krivov et al.ll 19981 ) and (|Mann et al.ll2000l ). the Lorentz 
force dominates the gravity for the dust grains smaller than 
0.1 ^.m i si z e, an d, as show their numerical simulations 
iMann et al.l (|200d ). the orbits of such grains get random- 
ized which is not th case for the grains with size larger than 
0.1 (im. Interestingly, for grains with the size smaller than 
0.01 Jim, the ratio charge-to-mass may be so high that they 
can be accelerated outward by the interaction with the so- 
lar wind as sho ws their detection by STEREO experiment 
iMever-Veriietet al.. (,2009i ). 
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A more detailed modelling of the presented data, taking 
into account the scattering geometry, dynamics, the density 
and size distributions of the dust grains, would be highly 
desirable to quantify further the dust physical properties 
near the Sun, but this is largely beyond the scope of the 
present article. 



5 CONCLUSIONS 

The reported here measurements of the Doppler shifts of the 
Fraunhofer lines, scattered by the dust grains in the solar 
F-corona, show that at the date of our observations the dust 
grains were on the orbit with a retrograde motion in a plane 
at i ~ 105°, i.e. nearly perpendicular to the ecliptics. This 
points to their cometary origin. Indeed, at the end of March, 
2006, SOHO recorded several sungrazing comets with the or- 
bital elements close to what was deduced from our measure- 
ments. We conclude that the contribution of comets to the 
dust content in the region close to the Sun can be important 
albeit variable in time. This contribution can explain the al- 
ready noticed change of the dust distribution from one in 
the axial symmetry far from th e Sun to tha t in th e central, 
or may be spherical, svmmetrv lMann et al.l ()2000l l. 

We also derive that the observed plane of the dust grains 
orbit is slightly different from that of the parent comet(s), 
which indicates that the size of grains is small, less than 
0.1 p.m, so that they are deviated from the initial orbit by 
the Lorentz force. This also means that the observed dust 
grains were released by the comet(s) shortly before our ob- 
servations. 

The importance of comets in a circumstellar environ- 
ment is general, let us recall e.g. the "Falling Evaporating 
Bodi es" recorded in the spectra of /3 Pic (e.g. iBeust et al.l 
Il998l ). It would be interesting to investigate whether they 
can provide a sufficient transport of the dust grains between 
far and close environments of the central star, and to con- 
tribute to the dissemination of crystalliz ed material in the 
recen tly detected exozodiacal dust discs l|Absil et al.ll200^ . 
I2009D . For a m ore detailed disc ussion on the possible role of 
comet see also lAugereaul (|2009l ) and references therein. 
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